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Abstract 


We derive a new set of sum rules between the neutrino mass and mixing 
parameters in vacuum and their effective counterparts in matter. The moduli 
of nine genuine lepton mixing matrix elements can then be calculated in terms 
of the matter-corrected ones, and the latter can directly be determined from 
a variety of long-baseline neutrino oscillations. We show that it is possible 
to reconstruct the leptonic unitarity triangles and CP violation in a similar 
parametrization-independent way. 
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I. INTRODUCTION 


Thanks to the Super-Kamiokande [1] and SNO [2] experiments, both solar and atmo¬ 
spheric neutrino oscillations are convincingly established. It turns out that neutrinos do have 
masses and lepton flavors are really mixed, just as expected in some grand unified theories. 
Two neutrino mass-squared differences and three lepton mixing angles have been measured 
or constrained by current neutrino oscillation data [1-5]. A more precise determination of 
these parameters has to rely on the new generation of accelerator neutrino experiments with 
very long baselines [6], in which leptonic CP violation may also be observed. The terrestrial 
matter effects in all long-baseline neutrino experiments must be taken into account, since 
they can unavoidably modify the genuine behaviors of neutrino oscillations in vacuum. 

To formulate the probabilities of neutrino oscillations in matter in the same form as those 
in vacuum, one may define the effective neutrino masses m* and the effective lepton flavor 
mixing matrix V in which the terrestrial matter effects are already included. In this common 
approach, it is necessary to find out the relationship between the fundamental quantities of 
neutrino mixing in vacuum (m, and V) and their effective counterparts in matter (m* and 
V). The exact formulas of m l and V as functions of m* and V have been achieved by a 
number of authors [7-10], and the similar expressions of rn l and V in terms of rhi and V 
have been derived in Ref. [11], The latter case is equivalently interesting in phenomenology, 
because our physical purpose is to determine the fundamental parameters of lepton flavor 
mixing from the effective ones, whose values can directly be measured from a variety of 
long-baseline neutrino oscillation experiments. 

This paper aims to reconstruct the genuine neutrino mixing matrix V and its unitarity 
triangles from possible long-baseline neutrino oscillations with terrestrial matter effects. Our 
work is remarkably different from the existing ones [7-11] in the following aspects: 

• We derive a new set of sum rules between ( rrii,V ) and frhi,V). It may perfectly 
complement the sum rules obtained in Ref. [12], We show that similar results can be 
achieved in the four-neutrino mixing scheme. 

• Our sum rules allow us to calculate the moduli of V ai (for a = e, /i,t and i — 1,2, 3) 
in terms of those of V a i. This approach, which results in some much simpler relations 
between \V Q i\ 2 and |I4i| 2 , proves to be more useful than that proposed in our previous 
work [10]. The sides of three unitarity triangles can also be derived in a similar way. 

• We express the neutrino oscillation probabilities in terms of \V a i\ 2 and the matter- 
corrected Jarlskog parameter J [13]. The fundamental lepton flavor mixing matrix 
V and its unitarity triangles can then be reconstructed from possible long-baseline 
neutrino oscillations straightforwardly and parametrization-independently. 

The only assumption to be made is a constant earth density profile. Such an assumption 
is rather reasonable and close to reality for most of the presently-proposed terrestrial long- 
baseline neutrino oscillation experiments [6], in which the neutrino beam is not expected to 
go through the earth’s core. 

The remaining parts of this paper are organized as follows. In section II, we figure out 
a new set of sum rules between ( rrii,V ) and ( rhi,V ) and discuss its extension in the four- 
neutrino mixing scheme. Section III is devoted to the calculation of \V a i \ 2 in terms of |I4,i| 2 . 
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We also show how to derive the sides of three unitarity triangles in a similar approach. In 
section IV, the neutrino oscillation probabilities are presented in terms of \V a i\ 2 and J . We 
illustrate the dependence of different oscillation terms on the neutrino beam energy and the 
baseline length. Section V is devoted to a brief summary of our main results. 


II. NEW SUM RULES BETWEEN V AND V 


In the basis where the flavor eigenstates of charged leptons are identified with their mass 
eigenstates, the lepton flavor mixing matrix V is defined to link the neutrino mass eigenstates 
(u \, v 2 , /Aj) to the neutrino flavor eigenstates (u e , v lt , v T )\ 


( v e 

'V 

\u T 


(v el 

Ve2 

V e3 \ 

( Vl \ 

K 1 

Vf/,2 

Ks 

v2 

VKi 

v t2 

v T 3 ) 

Us / 


( 1 ) 


A similar definition can be made for V, the effective counterpart of V in matter. The strength 
of CP violation in normal neutrino oscillations is measured by a rephasing-invariant quantity 
J (in vacuum) or J (in matter), the so-called Jarlskog parameter [13]: 
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( 2 ) 


where the Greek subscripts (a,/?, 7 ) and the Latin subscripts ( i,j,k ) run over (e,/i, r) and 
(1, 2, 3), respectively. 

The effective Hamiltonian responsible for the propagation of neutrinos in vacuum or in 
matter can be written as 

H - = h (M - Ml) = 21 ( VDlv ’) ■ 

7 = = _L {yblv t) , (3) 

where D u = Diag{m^, m|, and D u = Diag{mf, Wj}, E is the neutrino beam energy, 
M u and M v denote the corresponding genuine and matter-corrected neutrino mass matrices 
in the chosen flavor basis, m l and m* (for z = 1,2, 3) stand respectively for the neutrino 
masses in vacuum and those in matter. The deviation of H u from H u results non-trivially 
from the charged-current contribution to the u e e~ forward scattering [14], when neutrinos 
travel through a normal material medium like the earth: 

H v - H u = = Diag{a, 0, 0} , (4) 

Zrj 

where a = ^2GpN e with N e being the background density of electrons. Subsequently 
we assume a constant earth density profile (i.e., N e = constant), which is a very good 
approximation for most of the long-baseline neutrino experiments proposed at present. 
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Eq. (3) implies that {M v Ml) n = VDl n V t and ( M v Ml) n = VDl n V ] hold, in which 
n = 0, ±1, ±2, etc. To be explicit, we obtain 

(M„Mt) n aP = E (mfKaV^) , 

= E , (5) 

2=1 

where a and (3 run over e, // and r. The simplest connection between and ( M v Ml) n 

is their linear relation with n = 1; i.e., 

M v Ml = M V M\ + , (6) 

as indicated by Eqs. (3) and (4). With the help of Eqs. (5) and (6), one may easily find 

3 3 

E = E + akm . (?) 

2=1 2=1 

where A = 2 Ea. In the a f- /3 case, Eq. (7) reproduces the sum rules obtained in Ref. [12]. 
A new set of sum rules can be achieved, if the square relation 

[M v Mlf = ( M v M\f + + n v {M„Mt) + £% ( 8 ) 


is taken into account. Resolving Eq. (8) by use of Eqs. (4) and (5), we arrive at 

E (rhtV ai v£) = E {m 2 [m 2 + A (8 ae + M] } + A 2 8 ae 8 e/3 . 


2=1 


2=1 


( 9 ) 


Eqs. (7) and (9), together with the unitarity conditions 


E(t»b.) = E(K*V£) = 4s, 


( 10 ) 


2=1 


2=1 


constitute a full set of linear equations of V ai Vp t or V m Vp t for i — 1, 2 and 3. We shall make 
use of these equations to derive the concrete expressions of V ai V^ in terms of m iy fhi and 
VaiVpi in the next section. 

It is worth remarking that the sum rules in Eq. (9) are completely different from those 
presented in Refs. [10,15-17], where the relationship between 77,7 1 det. H v and 77“ 1 det H v 
has been used. Our new result can appreciably simplify the calculations of \V a i \ 2 in terms 
°f |Eil 2 (or vice versa), as one may see later. 

Eqs. (7) and (9) may be generalized to include the mixing between one sterile neutrino 
(u s ) and three active neutrinos (. v e , v jx and v T ). In this case, V is redefined to link the 
neutrino mass eigenstates (i/q, zq, v 2 , z/ 3 ) to the neutrino flavor eigenstates (zq, u e , u T ): 
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and V can be redefined in the same manner. The matrices D u , D u and in Eqs. (3) and 
(4) are now rewritten as 

D u = Diag {mg, m\, m 2 2 , m 2 } , 

D u = Diag {mg, m\, rh\, ^ 3 } , 

= Diag {A, A, 0,0} , (12) 

in which mo denotes the sterile neutrino’s mass, and A = \Z2GpN n E with N n being the 
background density of neutrons [18]. Different from A, A! measures the universal neutral- 
current interactions of u e , and v T with terrestrial matter. Both A and A' are assumed to 
be constant throughout this work. Then two sets of sum rules can respectively be derived 
from the linear and square relations given in Eqs. (6) and (8): 
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Of course, one may also make use of the relation between and (A4 u Ml) 3 to derive 

another set of sum rules between ( m^V) and (rhi,V), although the relevant calculations 
are somehow lengthy. It is then possible to establish a full set of linear equations of V a iVZ 
or VaiVpi f° r * = 0,1,2 and 3 in the four-neutrino mixing scheme, from which both the 
moduli of V a i and the sides of six unitarity quadrangles of V [19] can be derived in terms 
of A, A and the neutrino mass and mixing parameters in matter. Such an idea and its 
phenomenological consequences will be elaborated elsewhere. 


III. MODULI AND UNITARITY TRIANGLES OF V 


Now let us calculate the moduli of V ai by using Eqs. (7), (9) and (10) in the conventional 
three-neutrino mixing scheme, assuming that \ V a i \ 2 can directly be determined from possible 
long-baseline neutrino oscillation experiments. Those three equations are rewritten, in the 
a — (3 case, as follows: 
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(16) 


With the help of the relationship [15] 


E ™i = E m i + A > 

i =1 i=l 

we solve Eq. (14) and obtain the following exact results: 
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in which the neutrino mass-squared differences A 7 = m 2 — m 2 and A 7 = m\ — m 2 are 
defined. One can see that these results are more instructive and much simpler than those 
obtained in Ref. [10], because we have introduced a new set of sum rules in Eq. (9). 

Next we calculate V m Vg t in terms of E m E/ 3 * (for a. ^ (3). The former can form three 
unitarity triangles in the complex plane, which were originally named as A e with (a,/3) = 
(//, r), with (a, (3) = (r, e), and A r with (a, (3) = (e,/i) in Ref. [20]. Their effective 
counterparts in matter are then referred to as A e , A /t and A r . Taking account of Eqs. (7), 
(9) and (10), one may easily write down a full set of equations of E m V) 3 * with a ^ (3: 
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We solve Eq. (20) and arrive at 
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for A t , where A t] = mf—rrij denote the neutrino mass-squared differences in vacuum. These 
results are exactly the same as those obtained in Ref. [10], although a different approach 
has been followed. Eq. (22), (23) or (24) allows us to establish a direct relation between J 
and J defined in Eq. (2). A straightforward calculation yields JA 2i A 3 iA 32 = JA 2i A 3 iA 32 , 
which has been derived in Refs. [12,21] in a different way. 

For the sake of completeness, let us list the expressions of two independent A^ in terms 
of their effective counterparts in matter. 


A 31 - 
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where [11] 


x — A 21 + A 31 — A , 

V — A 2 iA 31 — A A 2 i (l — [1^21+ A 31 (l — |VA| 2 ) , 

1 2a : 3 — 9xy — 27AA 21 A 31 |V el | 2 

z = cos - arccos-^- . 

|_3 2 (x 2 - 3 yf 2 J 

Three independent A tJ can be expressed as 

An - ~ ^x 2 - 3 y z + yjz{l - z 2 ) , 

A 22 = - ^ 1 /x 2 - 3 y z - /3 (1 - z 2 ) - A 2 i , 

A 33 = + jj; zy/x 2 - 3 y - A 31 . 


(26) 


Note that A^ = A^ + Ajj = A tJ + A, /:i holds. Thus both A^ and A^ are fully calculable, 
once A 2 i, A 3i , A, \V e i\ and |V^ 2 | (or |V^ 3 |) are specihed. 

Note that the afore-obtained results are only valid for neutrinos propagating in vacuum 
and interacting with matter. As for antineutrinos, the corresponding results can simply 
be obtained through the replacements V ==>- V* and A =>- —A (and A' ==>• —A' for the 
four-neutrino mixing scheme). 


IV. LONG-BASELINE NEUTRINO OSCILLATIONS 

The matter-corrected moduli | VA-i | J and the effective Jarlskog parameter J can, at least 
in principle, be determined from a variety of long-baseline neutrino oscillation experiments. 
To be concrete, the survival probability of a neutrino v a and its conversion probability into 
another neutrino i/p are given by [15] 

P{y<x v a) = 1 — 4 (|ho*hA | 2 sin 2 Fji) , 

i<j 

P(y a - u p ) = -4 5: [Re (VciVpjV^VZ) sin 2 F jt ] - 8 J (sin F jt ) , (28) 

i<j i<j 

where (a,/3) run over (e, fi) x (n, r) or (r, e), Fji = 1.27A jiL/E with L being the baseline 
length (in unit of km) and E being the neutrino beam energy (in unit of GeV), and J has 
been defined in Eq. (2). The Re( V m Vpj VA V * i% ) term in P(u a —> Up) can be expressed as 

Re y*) = 1 (K k v; t \ 2 - \vj%f - fey;/) (29) 

with a ^ (3 and i 7 ^ j P k. Then we rewrite Eq. (28) as follows: 

RK - Va) = 1 - 4 (iVail&l 2 sin 2 F 21 + \V al V^\ 2 sin 2 F 31 + |V^ 2 XC 3 | 2 sin 2 F 32 ) , 

P{y a ^) = -2 [[VAiV^ii 2 ^! + \v a2 v; 2 \ 2 s 2 + |k 3 e; 3 | 2 -s 3 ] - 8JS, , 
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(30) 



where 


51 = sin 2 F 32 - sin 2 F 21 - sin 2 F 31 , 

5 2 = sin 2 F 3 i - sin 2 F 2 i - sin 2 F 32 , 

53 = sin 2 F 2 i - sin 2 F31 - sin 2 F 32 , 

S'j = sin F 21 sin F 31 sin F 32 . (31) 

Some comments on Eqs. (30) and (31) are in order. 

• Three oscillation terms in P{y a —> u a ) are associated with E m V)X| 2 (with i ^ j). 
Taking i = 1 and j = 2 for example, we obtain an effective unitarity triangle whose 
three sides are V e iV* 2 , V^iV* 2 and V T 1 V* 2 in the complex plane. This triangle was 
originally named as A 3 in Ref. [20]. Triangles Ai (for i = 2 and j = 3) and A 2 (for 
i = 1 and j = 3) can similarly be defined. Their genuine counterparts in vacuum are 
referred to as A* for * = 1,2 and 3. If the sides of Aj can all be determined from 
P(v a —> u Q ) in some disappearance neutrino oscillation experiments, it is then possible 
to calculate individual |VA| 2 and to extract V ol \ 2 by using Eqs. (17), (18) and (19). 

• Three sides of triangle A e , A M or A r are associated with three oscillation terms S) 
(for i = 1,2,3) of the appearance neutrino oscillation probability P{y a —> z/g), while 
the effective CP-violating parameter J is relevant to the oscillation term Sj. To 
illustrate, we plot the dependence of Si and Sj on the neutrino (or antineutrino) 
beam energy E in Fig. 1, where two typical neutrino baselines L = 730 km [6] and 
L = 2100 km [22] have been taken. The input parameters include A 2 i ~ 8 x 10” 5 eV 2 , 
A 32 ~ 2.3 x 10“ 3 eV 2 , $12 ~ 33°, d 23 ~ 45°, 0 V3 ~ 3° and <5 « 90° in the standard 
parametrization of V [15]. In addition, the terrestrial matter effects can approximately 
be described by A & 2.28 x 10” 4 eV 2 F/[GeV] [23]. Fig. 1 shows that S t may have 
quite different behaviors, if the baseline length is sufficiently large (e.g., L ~ 1000 km 
or larger). Hence the proper changes of E and (or) L would allow us to determine the 
coefficients of Si, i.e., \V ai Vpi\ 2 . The parameter J can also be extracted from a suitable 
long-baseline neutrino oscillation experiment, because the dependence of Sj on E and 
L is essentially different from that of Si. Provided all or most of such measurements 
are realistically done, the moduli of V m and leptonic unitarity triangles A e , A^ and 
A r may finally be reconstructed. 

• It is clear that both types of neutrino oscillation experiments (i.e., appearance and 
disappearance) are needed, in order to get more information on lepton flavor mixing 
and CP violation. They are actually complementary to each other in determining the 
moduli of nine matrix elements of V and its unitarity triangles. In practice, the full 
reconstruction of V from V requires highly precise and challenging measurements. A 
detailed analysis of the unitarity triangle reconstruction can be found from Ref. [24], 
in which the issues of experimental feasibility and difficulties have more or less been 
addressed. 

Let us remark that the strategy of this paper is to establish the model-independent 
relations between V and V, both their moduli and their unitarity triangles. Hence we have 
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concentrated on the generic formalism instead of the specific scenarios or numerical analyses. 
Our exact analytical results are expected to be a very useful addition to the phenomenology 
of lepton flavor mixing and neutrino oscillations. 


V. SUMMARY 

We have derived a new set of sum rules between (m*, V) and (m*, V) for both three- 
and four-neutrino mixing schemes. With the help of our sum rules, we have calculated the 
moduli of V ai in terms of those of V ai . The sides of lcptonic unitarity triangles A e , and 
A t can also be figured out in a similar way. We find that it is useful to express the neutrino 
oscillation probabilities in terms of \V Q i \ 2 and the matter-corrected Jarlskog parameter J. 
The fundamental lepton flavor mixing matrix V and its six unitarity triangles can then be 
reconstructed, at least in principle, from a variety of long-baseline neutrino oscillations in a 
straightforward and parametrization-independent way. 

Our analytically exact results will be applicable for the study of terrestrial long-baseline 
neutrino oscillations, from which much better understanding of the neutrino mass spectrum, 
lepton flavor mixing and CP violation can be achieved. 

One of us (Z.Z.X.) is grateful to CSSM at the University of Adelaide, where this paper 
was written, for warm hospitality. Our work was supported in part by the National Nature 
Science Foundation of China. 
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FIG. 1. Numerical illustration of four different oscillation terms appearing in P{v a —> vq) and 
P(V a —* Vp), where two typical neutrino beam baselines have been taken. 
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